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ABSTRACT 

We demonstrate that the South Pole kilometer-scale neutrino observatory 
IceCube can detect multi-TeV gamma rays continuously over a large fraction 
of the southern sky While not as sensitive as pointing atmospheric Cerenkov 
telescopes, IceCube can roughly match the sensitivity of Milagro. Also, IceCube 
is complementary to Milagro because it will observe, without interruption, a 
relatively poorly studied fraction of the southern sky. The information which 
IceCube must record to function as a gamma ray observatory is only the directions 
and possibly energies of down-going muons. 



1. Introduction 

One of the challenges of operating a high-energy neutrino telescope is to reject the 
background of cosmic ray muons traversing the detector. It is possible to use this background 
to do TeV-energy gamma ray astronomy. Showers generated by TeV gamma rays from 
sources in the southern hemisphere can produce muons that penetrate to the depth of the 
IceCube (Achterberg, et al. 2003) detector. Although relatively few muons are produced 
in a gamma ray shower, we will demonstrate that a directional source of TeV photons can 
produce a measurable excess over the background of muons produced by cosmic ray primaries 
at a rate of ~ 2KHz. By sheer statistics, a Crab-like object, the standard candle of TeV 
astronomy, will produce a yearly excess of 3-5cr, depending on the high-energy characteristics 
of the source. At TeV energies, sources burst to fluxes as high as ten times the flux of 
the Crab. We will show that IceCube's sensitivity to TeV gamma rays is comparable to 
Milagro, although neither can match the instantaneous sensitivity of the new generation of 
atmospheric Cerenkov telescopes. The latter, however, are only capable of observing a few 
degree patch of the sky on moonless, clear nights. With a South Pole location, IceCube is 
unique in that it observes the same sky without interruption, 24 hours per day, 7 days per 
week. It is most sensitive to sources overhead at the South Pole, a poorly studied portion of 
the southern sky. 
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While IceCube is a whole-sky muon and neutrino observatory at energies in excess 
of ~ 1 PeV, where the atmospheric background is negligible because of its steep energy 
spectrum, at TeV energies, the background of down-going cosmic ray muons prevents the 
detection of neutrinos from Southern sources. To operate as a gamma ray observatory, all 
IceCube must do is record the directions of the "background" muons and form a muon sky 
map of the Southern hemisphere. The muons will be triggered at a rate of close to 2KHz 
and reconstructed in real time. Eventually other shower properties, such as energy, could 
be recorded and cuts applied to identify muons of gamma ray origin in the much larger 
background of cosmic ray induced muons. 

TeV gamma ray astronomy has developed in several important ways in recent years. 
In the 1990's, galactic sources, such as the Crab Nebula, and extragalactic sources, such as 
Markarian 421, were observed to produce TeV photons. Today, 6 sources have been detected 
by multiple experiments with high significance. Another eight have been detected by a sin- 
gle experiment with high significance. While sources are few by astronomy standards, their 
high-energy behavior is fascinating and the accumulating data now defies a comprehensive 
theoretical explanation (Weekes 2003). After the pioneering work of EGRET (Mukherjee, 
etal. 1997 and Sreekumar, etal. 1997), Whipple (Petry, et al. 2002) and others, new in- 
strumentation is being developed and deployed such as Milagro (Atkins, et al. 2001) and 
the next generation of atmospheric Cerenkov telescopes (ACTs) including VERITAS (Ong, 
etal. 2003), MAGIC (Baixeras, etal. 2003), HESS (Hofmann, etal. 2002), CANGAROO 
(Kawachi et al. 2002 and Okumura et al. 2002) and others, as well as satellite based detec- 
tors such as GLAST (Kamae, et al. 2003). For a recent review of gamma ray astronomy, see 
Ong (Ong, 2003). 

Presently, only the Milagro experiment is capable of monitoring large portions of the 
sky for TeV-bright transient sources that may appear, such as gamma ray bursts and blazars. 
Milagro is, however, only capable of observing approximately a 90° diameter cone of the sky 
at one time, leaving many prospective transient sources unobserved in the TeV. 

2. Multi-TeV gamma rays In IceCube 

IceCube, now under construction, is designed to observe high-energy muons and show- 
ers generated in neutrino interactions inside of or nearby the detector volume. Multi-TeV 
showers created by photon primaries in the atmosphere will also generate muons, some of 
which may propagate to the depth of the detector. 

The number of muons with energy greater than produced in a shower with a photon 
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primary of energy E y is well approximated by 

N(E >E )~ 2 ' 14 X 10 ' 5 l - ^_ (1) 

Vl ^~ cosO {E^l 0086)^/0080)- [) 

6 is the zenith angle of the shower and all energies are given in units of TeV throughout 
this paper. This parameterization is based on the results of monte carlo and is accurate for 
muons with energies < E 7 /10 and > 0.1 TeV (Halzen, Stanev & Yodh 1996, Halzen, 
Hikasa & Stanev 1986). This parameterization describes muons from the decay of pions 
produced in the photoproduction of shower photons. Muons from muon pair production or 
the decay of charm particles, which are not important at the energies of interest, have been 
neglected. For a complete discussion, see (Drees, Halzen & Hikasa 1989). 

For a muon to be observed, it must reach the detector with energy above the threshold 
of the experiment. For IceCube, this threshold is ~ .1 TeV. As they travel from the surface 
to the detector, muons lose energy continuously according to 

^ = -a-PE, (2) 

where a = 2.0 x 10~ 6 TeVcm 2 /g and (3 = 4.2 x 10~ 6 cm 2 /g (Gaisser 1991 and Dutta et al. 
2000). The muon range is then 



a + f3Ef T 



(3) 



For a vertical and down-going muon to reach a depth of 1400 m with an energy of 0.1 TeV, 
the detector threshold, it must have an energy greater than ~ 0.56 TeV at the surface. This 
threshold is higher for sources not at the zenith as the required range for such muons is 
increased to R^j cos 9. Table 1 shows the surface threshold of muons as a function of zenith 
angle. 

The flux of a generic TeV gamma ray emitter can be parameterized as 

^ = ^10- 12 cm- 2 s- 1 . (4) 
dE^ E a v ; 

Again, E 1 is in units of TeV. Typical sources which we consider have values of a ~ 2. 

To calculate the number of muons observed from a given source, we convolve Eq. 1 with 
Eq. 4. The result is the number of muons with energy above the surface threshold from a 
source 

„ r /v n , Z^™ _ F 7 10 12 2.14 x 10" 5 E 1 

N,(> E,, sur ) ~ / dE,^- f— -, (5) 

JE^ min & ^Vsur (iV, SU r/COS0) 
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~ 2 x 10- 17 cm" 2 s^ 1 — - —In f , (6) 

cos#(E MiSur /cos#)« VE 7imin / 

where / is a correction factor, approximately given by 

Eq. 5 is only accurate when the integral is performed over many decades of energy Since 
this is not necessarily the case for us, / has been introduced to correct for this. This 
parameterization agrees with explicit monte carlo results (Halzen, Hikasa & Stanev 1986). 
We will use £< 7jm i n — 10 x i^ iSur neglecting the small probability that lower energy gamma 
showers produce detected muons. 

To determine if the number of muons calculated above is a significant observation, the 
background from cosmic ray muons must be understood. The background for IceCube is 
predicted to be approximately given by (Achterberg, et al. 2003) 

AW ^ 0.5(cosfl) 2 - 3 5 2 t (8) 

where 8 is the angular resolution of the detector. For IceCube, the angular resolution is 
expected to be better than 1°; we will take 5 = 1° as a conservative estimate. 



The signal to noise for a given source is then 

-In A cfr t / y/0.5 x 10 7 (cosfl) 2 - 3 5 2 t, (9) 



N sig 1.1 x 10- 16 F 7 (cos#) a - L53 , /E„ 



/1\T " pa-0.53 
V -'Vbkg -^sur 



1.6 x 10- 16 F 7 (cos#) a - 2 - 68 , /E. 



lnf^W^A (10) 



~ pa-0.53 

^/^sur s -^7,min 

where 5 is in degrees, t is the observation time in seconds and A cS is the effective area of the 
detector, 10 10 cm 2 in the case of IceCube. 

As a first example, we consider a hypothetical TeV-bright object similar to the Crab 
Nebula, but located at the south pole zenith. With an overhead source, the muon surface 
threshold is approximately 0.560 TeV. For the Crab, we use the values F 7 = 50.0 and a = 2.6. 
The signal to noise of such a source in IceCube is 

TV- / E \ 

Slg - 1.5 x In' 7 ' max ^ 



V 7 ^ 



bk g 



( 7 ' max y (in 

V 5.6 TeV/' V ; 



for one year of observation. If the gamma ray spectrum extends to 50 TeV, a 3.3-a observation 
is predicted. If the spectrum extends to 150 TeV, a 5-er observation is possible. This 
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calculation is somewhat conservative in our choice of -E 7 , m m = lO-E^sur- Including gamma 
rays below this energy would enhance the significance of the source. 

If a prospective source is not directly overhead, the significance of an observation is 
reduced. Considering again the Crab-like source, but with 9 = 20° and -E 7)max = 50 TeV, we 
find a significance of ~ 2.8-cr. For the same source, but with 9 = 45°, a 1.2-a is expected. 

There are several sources which are known to be, or are likely to be, TeV emitters in 
the southern hemisphere including the BL Lac object PKS 2155-304, supernova 1006, the 
pulsar 1706-44, the radio galaxy Cen-A, the supernova remnant RXJ1713-39, the starburst 
galaxy NGC 253 and the Vela Pulsar. There are many sources in the southern sky which 
have not been adequately studied at this time, in part due to the fact that the vast majority 
of TeV gamma ray experiments have primarily operated in the northern hemisphere. The 
CANGAROO (Kawachi et al. 2002 and Okumura et al. 2002) experiment has established 
TeV emission from a variety of sources in the southern sky at this time, but like any ACT, 
cannot possibly monitor large portions of the sky. For a list of hypothetical sources and 
IceCube's sensitivity, see table 2. 

In addition to continuous sources, Blazars such as Markarian 421 or Markarian 501 have 
been observed to flare to luminosities significantly greater than, for example, the Crab. If 
such an event were to occur in the southern sky, a flare of one week duration and ten times 
the flux of the Crab, could be observed at ~ 4.6-cr significance by IceCube. 

For particularly bright sources, it may be possible for smaller area neutrino detectors 
such as AMANDA (Andres, et al, 2001) or ANTARES (Aslanides, et al, 1999) to make 
such observations. Although these experiments have considerably less effective area, they 
also have significantly lower energy thresholds. For a summary of high-energy neutrino 
astronomy, see (Halzen & Hooper, 2002) or (Learned & Mannheim, 2000). 



3. Comparison With Milagro 

Unlike IceCube, Milagro is located on the surface and, therefore, has a considerably 
lower energy threshold. However, without the depth of an underground experiment, Milagro 
suffers from larger backgrounds. 

To establish the sensitivity of Milagro to TeV sources, we consider one of the sources 
which they have observed with high significance, the Crab Nebula. The gamma ray spectrum 
of the Crab is approximately 

<]Nl - 5 x 10- n £- 2 - 6 cm" 2 s- 1 TeV (12) 
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which corresponds to F y ~ 50 and a = 2.6 in our parameterization. 

We parameterize the effective area of Milagro as (Atkins, et al, 2001) 

( 4 x 10 6 E 1 - 39 cm 2 £ 7 > 1 TeV, 
A cff (£ 7 ) = I 4 x 10 6 E 2 - 35 cm 2 1 TeV > E 7 > .3 TeV, . (13) 
{ ~ cm 2 .3 TeV > E 7 

The number of signal events predicted from the Crab Nebula in Milagro is then given by 

/•-E 7 ,max ^jy 
^signal ^ / A^E^t— ^dE y (14) 

J e ahj^ 

~2xl0~ 4 xt / E} ( - 39 - 2M dE J + E 2 - 35 ~ 2 - 6 dE y ~ 6 x 10~ 4 events/sec (15) 

-JlTeV </.3TeV 

using -E 7jmax ~ 50 TeV. Over one year of observation, with the Crab observable in the sky 
approximately 20% of the time, we predict ~ 4000 events/yr in Milagro. This appears to be 
a slightly optimistic estimate, as Milagro has actually observed an excess from the Crab of 
approximately 3870 events/yr. 

The background rate observed in Milagro is approximately 2300 events/day for an object 
observable 20% of the time. After 1 year of observation with Milagro, we would expect a 
signal to noise from the Crab of 

Nsig 4000 

- ~4.4cr. (16) 



^/N^~ g V2300 x 365 

Using the rate actually observed by Milagro, 3870 events/yr for ~ 900 days, we get 

iVsig 3870 x 900/365 
— . ~ — =^- ~ 6.6 a, 17 

^/N^ V2300 x 900 

which is similar to the 6.4-cr result published by the Milagro collaboration (Atkins, et al., 
2003 and Atkins, et al., 2001). 

To emphasize the comparison with IceCube, for an overhead source with a ~ 2.5, 
Milagro is sensitive to objects with ~ 60% of the minimum flux required by IceCube. As 
the zenith angle increases, this factor becomes smaller, reaching ~ 20% by 6 = 45°. For a 
list of hypothetical sources and sensitivities for IceCube and Milagro, see table 2. 



4. Comparison With Atmospheric Cherenkov Telescopes (ACTs) 

Atmospheric Cherenkov Telescopes (ACTs) are designed to be highly sensitive to show- 
ers generated by TeV gamma rays. They must be pointed a source, however, and are not 
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capable of continuously observing large portions of the sky as experiments such as IceCube 
or Milagro can. 

For comparison, consider a generic ACT with a threshold energy of .25 TeV and effective 
area of 3.5 x 10 s cm 2 . The number of events observed for such a telescope is given by 



> 7 ,ma X dE 



r E -i- 

iVsig ^ A cS t (0.68) 2 / — '- F f 1(T 12 cm" 2 s" 1 (18) 

A cS t F 7 4.6 x 10- 13 
— 1 1 19 

F t 

~ 6.4 x 10" 4 - 7 - or r , (20) 

(a- 1) x 0.25 01 - 1 ' v ; 

where t is the observation time in seconds. The factor of (0.68) 2 appears due to the definition 
of angular acceptance in the estimation of the experiment's effective area. 

If the ACT has, for example, a 30% energy resolution and 0.001 steradians angular 
acceptance, the corresponding background consists largely of misidentified hadronic showers 
and electron-primary showers. This background can be approximated by 



AW ^ 3 x 10~ n 



+ 



0.08 



^7,th ^7,th 



A cS t (21) 



— > 0.13 events per second. (22) 
Combining these results, we get a signal to noise of 

N ■ F r- 

J?H. ~ 0.0018 x \ nr 1 x Vt. (23) 

To compare with an experiment such as IceCube and Milagro, consider again a hypothetical, 
Crab-like TeV-bright source with F 7 ~ 50 and a ~ 2.6. In less than 100 seconds of observa- 
tions time, a 5-cr observation could be made. Similarly impressive results are expected for a 
variety of sources. 

An experiment such as IceCube cannot possibly compete with ACTs in sensitivity to 
known continuous sources. Unlike ACTs, however, experiments such as IceCube or Milagro 
are capable of observing large portions of the sky and, therefore, are sensitive to transient 
sources such as gamma ray bursts or blazars. We will discuss such sources in the following 
section. 



- 8- 



5. Transient Sources 

In the previous section, we demonstrated that ACTs are far more sensitive to known con- 
tinuous sources than either IceCube or Milagro In this section, we will discuss the prospects 
for IceCube and Milagro to observe multi-TeV transient sources to which ACTs are not 
sensitive. 

Using the methods described in sections 2 and 3, we have calculated the signal to 
noise expected in both experiments for a variety of hypothetical transient sources. These 
results (see table 3) demonstrate Milagro's sensitivity to gamma ray bursts and blazars for 
values of F 7 of order 10 4 and possibly lower, depending on the spectral characteristics and 
duration of the burst. IceCube, while not as sensitive as Milagro in this respect, is capable of 
monitoring much of the, unobserved, southern sky for particularly TeV-bright transients. A 
short duration event (t ~ a few seconds), such as a gamma ray burst could be observed with 
5-o" confidence by IceCube for F 1 as low as 10 5 if its spectrum extends to very high energies 
and is located overhead. For a less ideal short duration source, only extending to tens of 
TeV, or at 9 ~ 30° — 40°, for example, could be observed only if F 7 ~ 10 6 , or so. Longer 
duration transients, such as blazars, which typically have flares of many minutes duration, 
can be observed with less fluence. 

Again, we would like to point out that no experiment other than IceCube will be capable 
of monitoring the southern sky for TeV-bright transients. 

6. Conclusions 

We have calculated the rates and sensitivities of IceCube as a TeV gamma ray observa- 
tory. IceCube, capable of detecting muons of ~ .ITeV and above, can observe the presence 
of muons generated in multi-TeV gamma ray showers, and distinguish these events from 
background given a sufficiently TeV-bright source. 

We have compared the sensitivities of IceCube, in this respect, to both Milagro and a 
generic ACT. ACTs are considerably more sensitive, but do not have the ability to monitor 
large portions of the sky continuously. Comparisons of IceCube with Milagro indicate that 
IceCube's sensitivity is comparable to Milagro's for sources near the southern zenith. Addi- 
tionally, these experiments are complementary, as they do not observe the same portions of 
the sky. If a TeV-bright transient source occurs in the southern hemisphere, IceCube may be 
the only experiment capable of monitoring it. No additional hardware or software is needed 
in IceCube beyond its planned design. 
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Table 1: The energy threshold for muons at the surface to reach a detector buried 1400 
meters below the Earth's surface in ice, as is the case for IceCube. Sources near the horizon 
have considerably high er thresholds for muons to successful ly propagate to the detector. 

Zenith Angle E^ SUT (TeV) 



0° 


0.561 


10° 


0.570 


20° 


0.601 


30° 


0.660 


40° 


0.765 


50° 


0.962 


60° 


1.391 


70° 


2.738 


80° 


16.546 



Table 2: The significance at which various sources could be observed by IceCube and MI- 
LAGRO in one year. All rates are calculated for sources with an a ~ 2.6 spectrum. For 
MILAG RO, we assumed that each source could be viewed ~ 20% of the tim e. 

Source Characteristics IceCube MILAGRO 





-- 1000,6 


no rp 


= 20 TeV 


38-o 


76-0 


F,= 


1000,0 


= 0° E 

' 7) max 


= 1000 TeV 


160-<7 


110-0 


F,= 


1000,0 = 


- 4^° F 


= 20 TeV 


12-0 


76-0 


F,= 


1000,0 = 


- F 

— *±<j , -u^max 


= 1000 TeV 


62-0 


110-0 


F,- 


= 100,0 : 


= 0° E 


= 20 TeV 


3.8-0 


7.6-o 


F ( = 


= 100,0 = 


= 0° B 

' J - J 7,max 


= 1000 TeV 


16-0 


11-0 


F,= 


= 100,0 = 


- F 

- *±<j , J-^^max 


= 20 TeV 


1.2-0- 


7.6-a- 


F,= 


100,0 = 


- F 

- *±<j , J-^'y ) max ~ 


= 1000 TeV 


6.2-o 


ll-o- 


F 7 


= 10,0 = 


- 0° E - 

' J - J 7,max 


: 20 TeV 


0.38-0 


0.76-0 


iV 


= 10,0 = 


--0° E = 

' J - J 7,max 


1000 TeV 


1.6-cr 


l.l-o- 


F ( - 


= 10,0 = 


AK° F 

*±<j , J-^ 7j max 


= 20 TeV 


0.12-a- 


0.76-0 


F,= 


= 10,0 = 


4^° F - 
*±<j , J-^'y ) max — 


= 1000 TeV 


0.62-0 


l.l-o- 



-13- 



Table 3: The observation time required to make a 5-cr observation in IceCube and MILAGRO. 
All rates are calculated for sources with ana~ 2.6 spectrum. 





Source Characteristics 


IceCube 


MILAGRO 


F,= 


= io 6 ,e = 


. no p 

- U , -C^max 


20 TeV 


0.55 sec 


0.027 sec 




io e ,e = 


0° E = 


1000 TeV 


0.031 sec 


0.012 sec 


F,= 


io 6 ,e = 




= 20 TeV 


5.5 sec 


0.027 sec 




io 6 ,# = 


P — 


1000 TeV 


0.21 sec 


0.012 sec 


F,= 


-- io 5 ,e = 




20 TeV 


55 sec 


2.7 sec 


F,= 


io 5 ,e = 


: 0° E = 


1000 TeV 


3.1 sec 


1.2 sec 


F,= 


w 5 ,e = 


/l 1 ^ P - 
*±tj , -i-J^max — 


= 20 TeV 


550 sec 


2.7 sec 


F,= 


io 5 ,e = 


/l 1 ^ P — 


1000 TeV 


21 sec 


1.2 sec 


F,= 


-- 10 4 , 6 = 


- n° p 

- U , -i-J^max 


20 TeV 


5500 sec 


270 sec 


F,= 


10 4 ,fl = 


w 5 J - J 7,max 


1000 TeV 


310 sec 


120 sec 


F,= 


10 4 ,# = 


4^° P 

*±tj , J-^7 l max 


= 20 TeV 


55, 000 sec 


270 sec 


F,= 


io 4 ,# = 


P — 
*±tj , J-^7 J max 


1000 TeV 


2100 sec 


120 sec 



